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£2 ' Abstract 

We calculate the initial density of relativistic global vortices (strings) produced 
at a quench and contrast it with the predictions of Kibble and Zurek. 
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Large-scale structure in the universe has been attributed |1[ to the production of 
cosmic strings (field vortices) in the phase transitions that are expected to occur in the 
Grand Unification era. However, it has been argued 0, |3|, || that there are similarities 
between the production of vortices in the early universe and in condensed matter systems 
( 4 He and 3 He in particular), where there is more accessible experimental data|| |^|. It 
is this suggestion that we shall be addressing in this letter, concentrating on one single 
property, the vortex density 

The simplest theory that displays vortices at a temperature quench is that of a 
complex scalar field = (0i + 20 2 )/\/2, determined by the 0(2)-invariant action 

s[</>] = J d*x(±(dM 2 - \m\t)4>l - ~A(0 2 J 2 )- (i) 

We identify the effective mass m(t) in Eq.p] with the equilibrium thermal mass, initially 
assumed positive. Once the quench is over, m 2 {t) takes its zero-temperature physical 
value — M 2 . In the interim, by reducing the temperature, we assume an explicit time- 
dependence m 2 (t) = —M 2 t/t Q , if —t Q < t < t Q and m 2 {t) = M 2 when t < —t Q (and is 
equal to — M 2 , for t > £q). Adopting different initial and final mass values (i.e. different 
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initial temperatures T = 0(M/y/X)) will have little effect, provided the initial mass is 
not too small. 

Estimates of the initial vortex density can be made in several ways. In Ref. || Kib- 
ble assumes that the equilibrium correlation length £ e q{t) = describes the corre- 
lation length £(t) of the field until time t ~ — t c , where £ e(3 (— t c ) = 1, on causal grounds. 
If £(£) is approximately frozen during the interval (—tc,t c ) then tc = tq 3 (2M)~ 2 ^ 3 , at 
which time the number of relativistic strings per unit area is 

n(tc) = 0(l/e eq (tc)) =o([At Q )-^ (2) 

where A = 0(M~ 2 ) is the cross-section of a classical vortex. The critical slowdown 
argument presented by Zurek|| sounds rather different, but the results are qualitatively 
the same in their power dependence on ig, supported by numerical simulations of the 
time-dependent Landau-Ginzsburg equations f?J, the standard dynamics of condensed 
matter systems. This early density is the largest as vortices subsequently annihilate each 
other but the details of this are irrelevant to our problem. The nonrelativistic counterpart 
of EqQ is in good agreement with the 3 He data||. Distinct vortices require n <C A^ 1 . 

In this letter we examine whether the result (0), based on the simple picture of 0, 
or the Landau-Ginzsburg equations of J7|, with its lack of direct reference to quantum 
fluctuations, is general. We shall see that the weakly-coupled quantum field theory of Eq.|l|, 
when used directly to estimate vortex densities, gives an initial density less than that of 
Eq.|2|, and qualitatively different. The reason is that the time tc is too early for a classical 
description of the field to be appropriate. This does not deny the usefulness of Eq.|| for 
condensed-matter systems. We first summarise our results and justify them later. 

It is because the formation of defects is an early-time occurrence that it is, in large 
part, amenable to analytic solution. Nonetheless, the production and freezing in of the 
vortices occurs over several time scales, as is most easily seen for an instantaneous quench 
(£q = 0), for which Eq.^Jis clearly inappropriate, a useful warning that the simple picture 
above has its limitations. However, once we have understood the instantaneous case the 
slow quench is not much more difficult. 

For the instantaneous quench there is, first, an interval Mt r = 0(1) before the field 
can respond. This is followed by an interval Mt sp = 0[ln(l/A)] during which domains 
form as a result of the growth of the long wavelength modes as they fall from the top of 
the potential hill towards its ground state value 0q = M 2 / A. After a further interval back- 
reaction (field self-interaction) will trap the field into one of its vacuum values, labelled 
by the phase a of ((f)), at each point in space. A consequence of this phase separation is 
that continuity and single valuedness will sometimes force the field to remain in the false 
ground-state at <fi = 0, the core of a tube of false vacuum, which characterises a vortex. 
By identifying string cores as lines of zeroes |§ of the field at the spinodal time t sp we 
shall find 

n ^ = °GtRI7A)) < 3 > 

and argue that this is, qualitatively, the density after the backreaction has taken effect. 

As the quench time £q is increased from zero we find different behaviour according 
as MtQ is 0(1), 0(ln(l/A)), 0(1/A) (or powers thereof), which we take to be qualitatively 
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distinct. Most simply, n of Eq.||is unchanged qualitatively as long as Mtq = 0(1). How- 
ever, once Mtq = 0[(ln(l/A)) Q ] with a > 1, or even a small power of 1/A, the situation 
begins to change. The response time t r now satisfies Mt r = 0((M£q) 1//3 ), and domains 
grow until a time given by Mt sp = 0[(Mt Q ) 1 / 3 (ln(l/A)) 2/3 ]. The corresponding den- 
sity is 



n = { (At Q mHi/x))^ ) <n{tc) - (4) 

Although this looks rather like n(tc) of Eq.||] appearances are deceptive, in that ln(l/A) 
is essentially a function of MtQ. 

We derive the results above by assuming Gaussian fluctuations. The vortex density 
n(t) at time t after the quench is started is given in terms of the equal-time Wightman 
function with the given thermal boundary conditions 

W/4|x-x'|;t) = (0 a (t,x)0 b (t,x')> = J fke ik - x P(k;t)5 ab . (5) 

P(k;t) is the power of the field fluctuations with wavevector k. Because of the non- 
equilibrium time evolution there is no time translation and, assuming no initial bias, W ab 
is 0(2) diagonal, W ab = S ab W. 

We count vortices by a judicious count of the line zeroes of the complex field. 
Care is needed since quantum fluctuations lead to line zeroes of the complex field on 
all distance scales. Let us cut off the field fluctuations at a wavenumber A. In our 
Gaussian approximation, with (<p a ) = at all times, this is equivalent to integrating over 
the 'environment' of short wavelength modes with k > A to construct a coarse-grained 
density matrix, from which W\ is derived. The resulting density of line zeroes is|| 

Of course, there will always be small-scale structure on vortices but only when the density 
is insensitive to A = 0(M), corresponding to their classical thickness, do we accept the 
line zeroes as representing vortices with a classical future. That is, we require the power 
k 2 P(k; t) in the fluctuations at A = 0(M) to be relatively small. 

We return to the instantaneous temperature quench (tg = 0) in which the effective 
mass m 2 (t) flips from +M 2 to — M 2 at t — 0. Initially, let us calculate the Gaussian 
W\{r \ t) for a free roll^, |10fl , built from the modes Xfc(^); 



£ + k> + m»(*) 



xi(t) = 0, (7) 



subject to xt(t) — e ±tuJint at t < 0, for incident frequency uji n = y/k 2 + M 2 . Imposing the 
KMS boundary conditions at t < for an initial Boltzmann distribution at temperature 
T determines W\(r;t) uniquely. For A = O(M), but A > M, say, W\(r;t) can be 
decomposed in an obvious way as W\(r; t) = W™(r) for t < 0, and 

W A (r; t) = Wt{r) + W w<M (r; t) + W A>w>M (r; t) (8) 
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for t > 0, where fTTH the second and third terms are the long and short wavelength fluctu- 
ations that grow from t = onwards. 

Initially, fluctuations are so large that the density depends critically on the scale 
A -1 at which we view the field. The equilibrium background term W™(r) has the form 

Wf(r)= [ fke^C(k), (9) 

J|k|<A 

where C{k) = coth(ui n (k) /2T Q ) /2uj in (k) encodes the initial conditions. The correlation 
length of the field is M~ l and, for cutoff A = O(M), the cold vortex thickness, we have 
riM = 0(M 2 ). There is a sea of interlocking line zeroes, separated only by the Compton 
wavelength. The oscillating field fluctuations make them unsuitable as serious candidates 
for string since, if we evaluate dn^/dA we find 



M dn A (t) 



n M (t) OA 



A=AI 



0(1). (10) 



A density of line zeroes so dependent on the scale at which it is measured cannot be 
compared usefully with that of (|2]). We term these the fluctuation vortices. 

For early positive times n A is equally sensitive to A but, once Mt 3> 1, the relevant 
term is W / |k|<M( r ; t), whose integral at time t is dominated by the nonperturbatively large 
peak|TT| in the power of the fluctuations at k around k , where tk\ ~ M. Once k <C M 



we have the required insensitivity. We find, approximately, that 

MT 

W W<M (r;t) ex p^e 2 ^exp{-r 2 M/4t}. (11) 

The prefactor comes from approximating C(fco) by Tq/M 2 . For t > the unstable modes 
grow until W|k|<M(0; t sp ) = O(0q) — 0{M 2 /X), determining the spinodal time t sp as 
Mt sp = 0(ln(l/A)). The long-distance behaviour of W™{r) = 0(e~ Mr ), with its correla- 
tion length £ = M~ l is a shortlived relic of the initial thermal conditions. After a time 
t r = 0(M~ l ) < t sp , it is rapidly supplanted by the behaviour of the expanding long wave- 
length modes W\^\ <M (r; t)/W\ k \ <M (0; t) w exp{-r 2 /£ 2 (t)} where £ 2 (t) = At/M w 4/fcg. 
With W\ k|<Af(0;t sp ) = 0(X x ) non-perturbatively large and W™(0) and W / A>|k|>Af(0; t sp ) 
of order A" 1 / 2 , 

n M (t sp ) « ~Jt—^ + 0(\ 1/2 ln(l/A))] (12) 

to give Eq.^|, as promised. The error term in the brackets is, in large part, a measure of 
the fluctuation vortices that we mentioned earlier. To check the stability of this density to 
changes in the coarse-graining scale A = O(M) a simple calculation shows that 9m {tap) — 
0{X 1 ^ 2 ln(l/A)) 1 for weak coupling. 

There are other, less direct, ways to understand why the strings only become well- 
defined once Wjw(0;i) is nonperturbatively large, even though the density is independent 
of its magnitude. Instead of a field basis, we can work in a particle basis and measure the 
particle production as the transition proceeds. Whether we expand with respect to the 
original Fock vacuum or with respect to the adiabatic vacuum state, the presence of a 
non-perturbatively large peak in k 2 P{k; t) signals a non-perturbatively large occupation 
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number iV& oc 1/A of particles at the same wavenumber fco[ll]. With Um of (|T^) of order 
k 2 . this shows that the long wavelength modes can now begin to be treated classically. 

^From a slightly different viewpoint, the Wigner functional only peaks about the 
classical phase-space trajectory once the power is non-perturbatively large [TO, |T3"[ from 
time t sp onwards. More crudely, the diagonal density matrix elements (field probabilities) 
are only then significantly non-zero for non-perturbatively large field configurations oc 
A -1 / 2 like vortices. 

So far we have done no more than estimate the density of coarse-grained line zeroes 
for a free roll at the time t sp at which nothing has frozen in. For this global 0(2) theory the 
damping of domain growth occurs by the self-interaction effectively forcing the negative 
m 2 (t) to vanish so as to produce Goldstone particles. It is instructive to adopt the simplest 
assumption that will lead to damping, in which we extend the time dependence of the 
free-roll mass to m 2 (t) = if t > t sp . The non-perturbatively large H^ic<m(^; t) then has a 
power spectrum that is the old one, frozen at time t sp , upon which oscillatory behaviour 
is superimposed. There is no qualitative change in the vortex density once the damping 
has taken place. 

To improve upon the free-roll result more honestly, but retain the Gaussian approx- 
imation, the best we can do is a mean-field approximation or, better, a large- N expansion 
(for O(N), N = 2). The main effect of the Goldstone production is to change the power 
in the field for very soft (k — > 0) modes, which has no significant affect on the position of 
the peak in the power spectrum, and hence the density flT|, |T2|1 . 

Whatever, but for varying coefficients of A, the modes Xk now satisfy the equation [p"Hj 



dt 2 



+ io 2 (t) + A / fpC(p)[xZ(t) Xp -(t) - 1] 



x* (*) = o. 



(13) 



where w 2 (t) = k 2 + m 2 (t). A rough estimate of the effects of the interactions prior to 
t sp can be obtained by retaining only the unstable modes in the integral for which the 
equations can be written 



d 2 l2 
df 2+k 



At) 



Xk(t) = o. 



(14) 



We adopt the approximate hybrid self-consistent equation for fj, 2 (t), 



H 2 (t) ~ M 2 - C\ 



(15) 



(C = 0(1)), which has the exponential growth of the WKB solution, but non-singular 
behaviour when fi(t) w 0. The exact solution to Eq.[15|for t < t sp is /i(t) = MtanhM(t sp — 
t). The condition fi(t sp ) = reproduces Mt sp = 0[ln(l/A)]. The first effect of the back- 
reaction is only to give a time-delay of 0(M _1 ), reducing the density by an unimportant 
0(1/M 2 ). 

For t > t sp oscillatory modes take over the correlation function to give oscilla- 
tions in W\(r;t) and its normalised derivatives and, through Eq.|| to oscillations in the 
string density. Interpolating the mode equation to larger t leads to oscillations with 
period 0(M~ 1 ). Explicit dissipative terms are not necessary to enforce damping since 
the O(((0 a (x) 2 )) 2 ) contribution to the mean-field Hamiltonian acts like a cosmological 
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constant [11]. In practice the backreaction rapidly forces /x 2 (t) towards zero if the cou- 



pling is not too smalllTT]]. Since the effect of imposing oscillations upon the frozen power 
spectrum is to leave the density, independent of normalisation, qualitatively unchanged 
in the nature of its A dependence it makes sense to identify the density at t sp with the 
density at freeze-in. 

All the results above were for the instantaneous quench. We now return to the 
original problem of slower quenches, in which the defects begin to freeze in by time 
t < tQ. For a free roll, with m 2 (t) as given before, the mode equation Eq.|^ is now of the 
form 

d 2 , 2 M 2 f 
+ k 2 



dt 2 tc 



X*(*)=0, (16) 



The exponentially growing modes that appear when t > t k = tqk 2 /M 2 lead to the 
approximate solution 

T / M \ 3 / 2 4Mt 3 / 2 



M\m(t)\ 

where £ 2 (£) = 2^/Uq/M. The provisional freeze- in time t sp is then, effectively, the time 
it takes to reach the transient groundstate 4>l(t) = —m 2 (t)/\. That is, W(0;t sp ) = 
0(<P 2 (t sp )), giving, for Mtq < (1/A), 

Mt sp ~ (Mt Q ) 1 / 3 (ln(l/A)) 2 / 3 ~ Mt c (ln(l/A)) 2 / 3 , (18) 

greater than Mtc- 

At this qualitative level the correlation length at the spinodal time is 

M 2 e{t sp ) * (Mt Q ) 2 / 3 (ln(l/A)) 1 / 3 . (19) 

The effect of the other modes is larger than for the instantaneous quench, giving, at t = t sp 

n M = -^[1 + 0(A 1 / 2 (Mt Q ) 4 / 3 (ln(l/A))- 1 / 3 )]. (20) 

The error term is due to fluctuation vortices. It is small if Mtq is a very small power 
of A -1 . However, it is more comfortable to take Mtg = 0((ln(l/A)) a ) (a > 1) or, 
equivalently, as a power (greater than one) of Mt sp , still permitting Mtq > Mt sp . 
Then for small enough A, we recover tim of Eq.^. As could be anticipated from Eq.|20|, 
©m(*s P ) ~ 0(A 1 / 2 (Mtg) 4 / 3 (ln(l/A)) -1 / 3 ), small enough to give an identifiable network 
on freezing in provided MtQ is constrained as above. Equivalently, by time t sp we have 
nonperturbatively large occupation numbers for the long wavelength modes, as evidenced 
by the WKB solutions. 

The density of Eq.§ is larger than that of Eq.^0] by the factor (ln(l/A)) 1/ ' 3 which 
is, in principle, a large number even if in practice it is only a small multiple. However, 
EqfJ is a misleading upper bound when considering definition of the network since a 
simple calculations shows that, at time t c , ©Af(^c) = 0(1). At this time, with no large 
occupation numbers, a classical description of strings is inappropriate. The bound set by 
causality is not saturated in our weak-coupling theory. 
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For t < t sp the backreaction is more complicated than for the instantaneous quench. 
A rough approximation is obtained by iterating the free- field behaviour. The mode Xk{t) 
then grows until time t^, defined by 



MHt CXT ( M \ 3 / 2 

The spinodal time t sp is the solution to Eq.|21] for k = 0. In terms of t sp the mode 
Xkif) stops growing at a time At(k) before it would have stopped had the backreaction 
been instantaneous, where MAt(k) = 0((k 2 /M 2 )(tQ/t sp ) 3 ^ 2 ). From Eq.|TT| this can be 
reinterpreted as a timelag At, where MAt = (^q/^sp) 1 ^ 2 - Although this is logarithmically 
large, relatively to t sp it is the same as for the instantaneous quench. We have to treat 
these ratios with care, since there is no self-consistency in them, but they continue to 
suggest that the initial effect of the back-reaction is benign. 

It seems that the initial effect of the backreaction after time t sp is to set up oscil- 
lations of period 0{M~ l tq/t sp ) followed, once fi 2 (t) « 0, by further damped oscillations 
over a time scale St = 0((^/tQt sp / M) l l 2 ). As for the instantaneous quench, the density 
at t sp determines the density at freeze-in. 

The uncertainty in string length because of the fluctuations in the short wavelength 
modes persists at the same level, increasing as the quench rate is slowed. Once the rate is 
slow enough the WKB approximation manifestly breaks down. This is almost certainly 
the case if the growing modes catch up with the moving minima within the Ginzburg 
regime by time to = O(Xtq). Simple inspection shows that fluctuation vortices dominate 
the density and m (to) is large. With no strong peaking in the power k 2 P{k) and no 
large occupation numbers that betoken classical behaviour the strings do not provide a 
well-defined network. This is despite their low density, which can be made compatible 
with the equilibrium fluctuation density uq = 0(AA _1 ) that was originally considered to 
be the relevant density for defect formation [|l|, PJ] . 
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